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Super Efimov effect of resonantly interacting fermions in two dimensions 
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We study a system of spinless fermions in two dimensions with a short-range interaction fine- 
tuned to a p-wave resonance. We show that three such fermions form an infinite tower of bound 
states of orbital angular momentum I = ±1 and their binding energies obey a universal doubly 



exponential scaling, £3 oc exp(— 2e 



37rn/4+9\ 



at large n. This "super Efimov effect" is found 



by a renormalization group analysis and confirmed by solving the bound state problem. We also 
provide an indication that there are I = ±2 four-body resonances associated with every three-body 
bound state at E^ oc exp( — 2e 37rn,/4+8_0 ' 1S8 ) . These universal few-body states can be observed in 
ultracold atom experiments and should be taken into account in future many-body studies of the 
system. 
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Introduction 

Recently topological superconductors have attracted 
great interest across many subfields in physics 0, 
This is partially because vortices in topological super- 
conductors bind zero-energy Majorana fermions and obey 
non-Abelian statistics, which can be of potential use for 
fault-tolerance topological quantum computation 0, Q . 
A canonical example of such topological superconductors 
is a p-wave paired state of spinless fermions in two dimen- 
sions which is believed to be realized in Sr 2 Ru0 4 0. 
Previous mean-field studies revealed that a topological 
quantum phase transition takes place across a p-wave 
Feshbach resonance 0-Q ■ 

In this Letter, we study few-body physics of spinless 
fermions in two dimensions right at the p-wave reso- 
nance. We predict that three such fermions form an infi- 
nite tower of bound states of orbital angular momentum 
£ = ±1 and their binding energies obey a universal dou- 
bly exponential scaling, 



E. 



(«) / 
3 oc exp( 



(1) 



at large n. Here 9 is a non-universal constant defined 
modulo 37r/4. This novel phenomenon shall be termed 
a "super Efimov effect" because it resembles the EA- 



TABLE I. Comparison between the Efimov effect versus the 
super Efimov effect. 



Efimov effect 



super Efimov effect 



three bosons 
three dimensions 
s-wave resonance 

e = o 

exponential scaling 



three fermions 
two dimensions 
p-wave resonance 
t = ±1 
doubly exponential scaling 



mov effect in which three spinless bosons in three di- 
mensions right at an s-wave resonance form an infinite 
tower of £ = bound states whose binding energies obey 
the universal exponential scaling, E^ 1 ' oc e ~ 27r "/ s o ^ w ith 
so ~ 1.00624 [lCj] (see Table Q] for comparison). We also 
provide an indication that there are I — ±2 four-body 
resonances associated with every three-body bound state 
at 



Ef ] ocexp(-2e 37 ™/ 4+e -°- 188 ), 



(2) 



which also resembles t he p air of four-body resonances in 
the usual Efimov effect Hlll2j- These universal few-body 
states of resonantly interacting fermions in two dimen- 
sions should be taken into account in future many-body 
studies beyond the mean-field approximation. 



Renormalization group analysis 

The above predictions can be derived most conve- 
niently by a renormalization group (RG) analysis. The 
most general Lagrangian density that includes up to 
marginal couplings consistent with rotation and parity 
symmetries is 

C = V f (idt + V' + ^ (idt + ^- - eo 
+ g4>li}{-iX7 a )vb + g4>\-iX7^ a )^(i) a 

+ V 3 1^<t>\M + V4 4>\4L a (t>-a<t>a + v'i fiifiifiafia- (3) 

Here and below, % = m = 1, V± = ± iS/ y , and sums 
over repeated indices a = ± are assumed, ip and <f>± 
fields correspond to spinless fermion and I = ±1 compos- 
ite boson, respectively. The p-wave resonance is defined 
by the divergence of the two-fermion scattering ampli- 
tude at zero energy, which is achieved by tuning the bare 



FIG. 1. Self-energy diagram of (j> a field. Solid and dashed 
lines represent propagators of ifr and <f> a fields, respectively. 



FIG. 3. Diagrams to renormalize the four-body couplings v 4 
and v' 4 . 



detuning parameter at £o = g A /(2tt) with A being a 
momentum cutoff. 

Since the other couplings (<?, U3, U4, v'4) are all dimen- 
sionless, our effective field theory (|3]) is renormalizable 
and its renormalization can be performed in a similar 
way to Ref. [l3j]. The self-energy diagram of (j) a field 
depicted in Fig. Q] is logarithmically divergent and renor- 
malizes the two-body coupling g. The RG equation that 
governs the running of g at a momentum scale e _s A is 
found to be 



which is solved by 



ds 



' 2tt'' 



(4) 



(5) 



ff 2 (0) 



Accordingly, the two-body coupling is marginally irrele- 
vant, i.e., it gets weak toward the infrared limit s — > 00. 

We now turn to the renormalization of the three-body 
coupling W3. Diagrams that renormalize V3 are depicted 
in Fig. [2] By taking into account the contribution from 
the 4> a field renormalization, the RG equation that gov- 
erns the running of V3 is found to be 



dv3 
ds 



I6g llg 2 V3 



3tt 



3tt 



24 
3tt ' 



(6) 



In the infrared limit s — > 00 where the two-body coupling 
scales as g 2 — > tt/s [141 ] . the analytic solution to Eq. © 
becomes 



v 3 (s) -> 



2tt 



1 



col; I -(In s 



(7) 



where the angle 9 depends on the initial conditions of g 
and ^3 at the ultraviolet scale s ~ 0. 





FIG. 2. Diagrams to renormalize the three-body coupling v$. 



We find that SU3 at large s is a periodic function of 
Ins and diverges at Ins = 37rn/4 + 9. These divergences 
in the three-body coupling indicate the existence of an 
infinite tower of energy scales in the three-body system 
in the £ = ±1 channels. Following the usual Efimov 
effect [TH, we identify these scales with binding energies 
which leads to the prediction presented in Eq. (JTJ). This 
identification will be confirmed later by solving the bound 
state problem. 

Similarly, the four-body couplings i>4 and v\ are renor- 
malized by diagrams depicted in Fig. [3] By taking into 
account the contribution from the <p a field renormaliza- 
tion, the RG equations that govern the running of V4 and 
v\ are found to be 



(8a) 



dv4 8g 4 


2g 2 v 3 


2g 2 v 4 




ds TT 


TT 


TT 


TT 


dv' 4 4 5 4 


2.9 2 f 3 


2g 2 v> 4 




ds TT 


TT 


TT 


TT 



(8b) 



respectively. Numerical solutions to these RG equations 
are shown in Fig. |4] and both SV4 and si>4 at large s are 
periodic functions of Ins with period 37r/4. We find that 
t>4 corresponding to the t — channel diverges only at the 




FIG. 4. RG evolutions of SV4, sv'4 (solid curves) and SV3 (dot- 
ted curve) with initial conditions g(0) = 10, 1)3(0) = ilO -6 , 
v 4 {0) = v' 4 (0) = 



FIG. 5. Two-fermion scattering with the amplitude repre- 
sented by the double line. 



FIG. 6. Three-fermion scattering with the amplitude repre- 
sented by the bulb. 



same points as 1)3. On the other hand, v\ corresponding 
to the I = ±2 channels has an additional divergence at 
Ins = 37m/4+0— 0.188 associated with every point where 
«3 diverges. Following the RG study of spinless bosons in 
three dimensions fl6j j , we identify these additional diver- 
gences with four-body resonances. We checked that this 
result is independent of the choice of initial conditions 
which leads to the prediction presented in Eq. 



Bound state problem 

The above predictions derived by the RG analysis can 
be confirmed by solving bound state problems. We em- 
ploy a model Hamiltonian 



f dk k ,+ v-^ f dkdpdq 

X Xa(p)x-a(q)^{ +p ip{_ p ipk^ q iPk +q , (9) 



which describes spinless fermions interacting by a sepa- 
rable p-wave potential. The form factor 



(10) 



is introduced to regularize the ultraviolet behavior. By 
summing diagrams depicted in Fig. [5j the two-fermion 
scattering T-matrix with incoming q and outgoing p rel- 
ative momenta is obtained as 



T(E;p,q) = 



167r|p||q| cos(tp p - ip q ) e 



-(pW)/(2A 2 ) 



fs - A 2 -Ee- E / A2 Eii-E/A 2 ) 



(11) 



where E = lu — k 2 /4 + i0 + is the total energy in the 
center-of-mass frame and Ei(w) = J dte^/t is the ex- 
ponential integral. Accordingly, the p-wave resonance 
is achieved by tuning the bare two-body coupling at 
v Q = 2tt/A 2 . 

A three-fermion scattering problem can be solved in 
a similar way to the corresponding problem in three di- 
mensions [IH [l9| • The T-matrix elements satisfy coupled 
integral equations depicted in Fig. [6] When the total en- 
ergy approaches a binding energy E —n 2 < 0, the T- 
matrix factorizes as T a b(E;p, q) — > Z a (p)Z£(q)/(E + K 2 ) 



and the resulting residue function satisfies 



Za(p) 



dq (p + 2 g )_ a e-( 5 P 2+5 <J 2 + 8 P-J)/(8A 2 
2^ 



P 



pq 



E b =±(2P + <l)bZ b (q) 



(Is 2 



+ K 



K 2 )/A 2 ) 

(12) 



It can be seen that Z + (p) = e 1 ^ 1 ' tPp z + (p) couples with 
Z-(p) = eW+Vvrz-ip). We focus on the I = +1 channel 
in which the super Efimov effect is expected to emerge. 

We first solve the two coupled integral equations 
(|12p analytically within the leading-logarithm approxi- 
mation [20|. We assume that the integral is dominated 
by the region where k <C q <C A and split the integral into 
two parts, k <C q <C p and p < 5 < A, in which a sum 
of p and q in the integrand is replaced with whichever is 
larger. This approximation simplifies Eq. (|12[) to 



; z +(p) 



f«-(p) 



p dq z+(q) j" A dq z+(q) + Z-(q) 
q \nk/q J p q 



In A/q 



p dq z+(q) 
q In A/9' 



(13a) 
(13b) 



By changing variables to x = lnlnA/p and y = In In A/q 
and defining A = lnlnA/re and C±( x ) — z ±(p)i we obtain 



C-(s) 



dyC+(y)- I dy[( + (y) + C-(y)} 

(14a) 



dy(+{y)- 



(14b) 



Then the differentiation of Eqs. (TP4"]) with respect to x 
results in two coupled differential equations 

3 



-C' + {x) = -C-(x) 



-C(x) = C + (x) 



(15a) 
(15b) 



with boundary conditions 
3 



C+(-oc) = - / dy(+(y) 



2 C-(A)=0. 



(16a) 
(16b) 



4 



TABLE II. Lowest five binding energies Eg = —Kn/m in the 
form of \ n = lnlnA/«; n . 



n 


Yi 


A n — A„_i 





0.5632 




1 


2.770 


2.207 


2 


5.078 


2.308 


3 


7.430 


2.352 


4 


9.785 


2.355 


oc 




2.35619 



The differential equations ([15]) with the boundary condi- 
tion (|16b[) are solved by 



C+O) = cos( -(a;- A) 



C-(x) = sin( -(a; - A) ) . 



(17a) 
(17b) 



The other boundary condition (|16al) constrains the al- 
lowed value of A, while we cannot determine its value 
within the present approximation because it is sensitive 
to the ultraviolet physics at q ~ A. However, owing to 
the periodicity of solutions (fTT)) . if A = 9 is a solution, 
then A = 3-7m/4 + 9 must be all solutions, which is con- 
sistent with our previous RG analysis. We note that the 
solutions for £ = — 1 are obtained simply by exchanging 
+ •<->•—, while the same approximation does not yield 
any solution for £ ^ ±1. The double- logarithm scaling of 
solutions was also found in studying a scattering problem 
of three fermions [jlj . 

We also solved the two coupled integral equations (fT2|) 
with £ = ±1 numerically and obtained binding energies 
of three fermions shown in Table HU We find that they 
obey the universal doubly exponential scaling with period 
37r/4 « 2.35619, which indeed confirms our prediction in 
Eq. While our prediction for £ = ±2 four-fermion 
resonances in Eq. ([2J can be confirmed in principle in a 
similar way to the corresponding problem in the usual 
Efimov effect 



study. 



Ill Il2| , we defer this problem to a future 



Conclusion 



In this Letter, we predicted the super Efimov effect in 
which resonantly interacting fermions in two dimensions 
form an infinite tower of I = ±1 three-body bound states 
whose binding energies obey the universal doubly expo- 
nential scaling as in Eq. ([1]). We also argued that there 
are £ — ±2 four-body resonances associated with every 
three-body bound state as in Eq. ©. It is possible to 
extend our analysis away from the p-wave resonance and 



draw a binding energy diagram as a function of the de- 
tuning in analogy to the usual Efimov effect [22j. These 
universal few-body states should be taken into account 
in future many-body studies of the system such as the 
topological quantum phase transition [jj @] ■ Our super 
Efimov states can be observed in ultracold atom experi- 
ments with p-wave Feshbach resonances [23| in the same 
way as the usual Efimov states through a measurement of 
three-body or four-body recombination atom loss 24 , 2|| . 
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